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Employing  the  Karmen  Method  for  Calculating  a  Turbulent 
Boundary  layer  on  a  Plata  in  a  Gaa  Flow 

*7 

L.  0*  Loytsanskiy;  Tu#  V.  Lapin 
Introduction 

In  the  theory  af  turbulent  boundary  layer  in  a  compressible  fluid  at  greater  velo 

cities  two  historical  trends  have  been  noticed*  The  first  one  of  these9  discovered  by 

semiemplrlcal  / 

F.I*Frankl  and  V*V*VoyBhel[iJ  was  based  on  direct  transfer  of  the/  theo- 

ry  by  Karmen  into  gas  dynamic  formulae*  Van  Oraat[2^  followed  an  analogous  path9but 
with  the  application  of  Rrandtl  formulae# 

In  reports  of  the  second  trend  were  used  known  Dorodnitayn^J  variables  and  it  has 
been  assumed,  that  the  corresponding  fornulae  of  the  eemiempir ical  theory  of  turbulences 
should  be  compiled  in  these  variables  # 

The  answer  to  the  question  on  which  trend  offers  results  more  closer  to  the  truth9 
could  be  given  only  by  experiments  at  greeter  M-numbersv  but  by  the  time  of  sppearanoe 
of  both  trends  experimental  data  were  still  lacking*  At  M-numbera,  not  much  exceeding 
one9  both  trends  have  given9  naturally,  results  slightly  differing fr on  each  other* 

Only  relatively  recently,  after  data  have  been  obtained  on  the  friction  resistance  of 
flat  surfaces  at  greater  M- numbers  (reaching  up  to  6  ,  it  became  clear,  that  only  the  first 
trend  offers  proper  results# 

W#  would  like  to  point  out,  that  in  the  Frank  1/Voyshel  report  the  authors  stood 
along  the  road  of  directly  generalising  the  Kerman  method,  having  simplified  only 
with  the  assumption  of  friction  stress  constancy  eeroee  the  boundary  layer*  Following 
this  path,  they  at  first  find  the  fora  of  the  profile  of  velooity  in  the  oross  eeetioa 
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of  the  layer *  then  in  an  ordinary  way  they  find  the  a o-called  law  of  "Reels, 

tivity"*  i.e.relatiooship  between  local  coefficient  of  friction  and  R-aumber  of  the 
boundary  layer  (R-number  *  Reynolds  or  Re-number  )#  Then*  excluding  this  Re -number 
from  equation  of  the  "law  of  resistivity  •  and  from  equation  of  pulses,  they  obtain 
the  sought  for  relationship  between  the  coefficient  of  resistance  and  the  Re-number* 
formulated  by  the  velocity  of  an  oncoming  flow  and  the  abscissa  of  given  point  on  the 
plate* 

The  degree  of  approximation  9  adopted  by  Frankl  and  Voyshel,  enabled  them  to  per¬ 
sona  lly  carry  out  calculations  but  only  up  to  I-t- numbers, slightly  exceeding  one*  The 
change  over  to  greater  M-numbcra  requires,  apparently,  either  an  even  greater  complexity 
and  besides  that  a  complex  calculating  method  by  these  authors,  or  direct  application 
of  numerical  methods* 

It  is  shown  in  this  report  that  when  applying  the  Karmen  method  and  assuming  frie~ 
tion  stress  constancy  across  the  boundary  layer*  there  is  a  simpler  way  of  formulating 
a  solution  ^requiring  no  introduction  of  the  concept  concerning  the  Reynolds  number  of  the 
boundary  layer  and  compilation  of  the  "law  of  resistivity  ",  This  method  considerably 
simplifies  the  investigation  of  problems  pertaining  to  a  turbulent  boundary  layer  in  a 
gas  flow*  The  use  of  a  simple  asymptotic  break  down  allows  comparatively  easily  to 
generalize  the  Karmen  theory  of  plate  resistance  in  an  incompressible  liquid  for  the 
case  of  a  gas  flow  with  greater  M-numbera. 

1*  Turbulent  Friction  on  a  Plate  in  a  Flow  of  Incompressible 
Liquid* 

The  integral  equation  of  pulses  for  a  pi.  te, around  which  an  incompressible  liquid 
flows  at  zero  angle  of  attack*  has  the  form  of  t 


where  z  -coordinate ,directed  along  the  surface  of  the  plate;  and  -  velocity 

and  density  of  the  oncoming  flow;  -  friction  stress  on  the  vall;^**-  thickness  of  pulse 

loss*  expressible  by  eqoalltys 
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In  expression  (l#2)  we  substitute  the  unlTersal  coordinates! 

’  t-  [1.3) 

here. .*.7  •coordinate  along  the  noraal  to  the  surface  of  the  plate  1  *  coefficient 

of  kinematic  viscosity  $  rm  »  ~l>  ~Tc7  •dynamic  Telocity. 

1  Q  - 

Then 

»•'  Q  t) 

0 

where 

(1.5) 

Equality  (l*4)  can  also  be  written  in  fora  of  s 

-IzH-l'-iMi)-  [l.t) 

u  ' 

ct>i  . 

The  value  ^  *  ~<f*f  1-11  ths  ^UBinar3r  sublayer  is  equal  tc  one  (^sl).  In  the 

turbulent  nucleus  to  determine  0?  we  use  the  Karman  ratiet 

where  x  -  turbulence  constantj^qualling  0*4  (the  hatchure  sign  in  the  equality  (1*7) 
designates  a  derivative  by  y)#  Adopting  a  simplified  assumption  about  the  friction 
stress  constancy  across  layer  ^  m*X&  and  transition  in  (l«7)  toward  universal  coot 
dinates  Cl«3)  will  obtains 

it*) 

The  minus  sign  appears  after  ext  action  of  the  root  in  connection  with  the  feet, 
that  y*  In  equation  (1*8)  we  change  the  argum&t  and  function  places,  then  we 
obtain  an  equations 


V  -  (1.1) 


from  which  it  is  svidsnt  that 


(l-lo) 


dtf 


It  vs  should  designate  by  the  lsttsr  f  ths  aagaituds  of  the  dsrlvatlvs  ou 
ths  boundary  of  ths  lamiuary  sublaysr  a  from  ths  sids  of  ths  turbulsnt  nuelo 

us,  and  writing 
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then  • qua t Ion  (1*10)  will  aequlra  the  fora  of 

.  {1/2) 

whore  a  -  turbulence  c  onsta  at  ^equalling  11*5# 

Returning  to  equation  (1.6)  we  want  to  point  out,  that,  strictly  speaking,  the 
integral  in  the  right  part  of  this  equation  should  have  been  broken  down  inte  two 
sections*  0  a  and  a  h  and  place  in  teach  one  of  these  a  value  .  naiaely 

^  *  1  in  the  first  and  according  to  Equation  (1.12)  •  in  the  second.  But  in  view  of 
the  relatively  thin  (relative  thinness)  of  the  leminary  sublayer  it  is  possible  to 
bypass  the  first  section,  extending  the  second  one  to  the  vail.  As  is  shown  by  calcn-* 
lations  in  first  approximation,  at  greater  h  no  noticeable  difference  is  obtained. 

In  this  way,  eliminating^  fran  equation  (l.6)  through  (1.12)  we  will  obtain  the 
following  expression  for  pulse  loss  thickness t 

6**  - 


--)**£.  (//  3j 


-  <¥ 
where  u  *  - 


The  value  xh,  as  la  evident  from  its  determination,  considerably  exceeds  one.  In 
these  cases  it  is  convenient  to  uee  the  following  general  idea  of  the  integral  (con¬ 
taining  the  exponential  function  ),  which  is  obtained  as  result  of  integrating  by 


parts i 


Kh 


/«>  --iS- 


JjV> 


ait) 


*n  thia  otM(  vbea  tba  aarita  in  aquara  parenthaaaa  ramalna  undiarupted,  fonoula 
(1.14)  can  ba  eonaidarad  aa  an  aaymptotie  aarioa,  axpraaaing  the  ralua  of  the  integral 
at  graatar  h« 

Making.  (1*13)  an  Integra t ion  with  tba  aid  of  (1*14)  va  will  obtain  an  axpreaaio. 
for  polaa  loaa  thlakaaaa  t 

V#  will  also  obtain  an  axpraaalc.  for  tba  fon>»paraJMtor  H| 
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u  6* 

H  -  -6T 


(//6) 


where  <)  +  •  thickness  of  squeezing  out  the  boundary  layer*  expressed  by  equation 


»•  -J'c (/  I'l) 


Using  a  method  ^analogous  to  the  one  shown  above*  it  is  possible  to  find  the  fol~ 
lowing  expression  for  the  thickness  of  squeeze  out < 


6*  =,±i^Le*h. 


{1.18) 


Compiling  the  ratio  according  to  equations  (1*1$)  and  (1*15)  v*  *111  obtain 

a  tens  for  Hs 


-r  (//*; 


which  at  Re  numbers  10  -  10*  gives  a  value  H  =  1*3*  well  conforming  with  the  experiment* 
Turning  to  expression  (l*15)t  ve  went  to  mai  tion*  that  in  view  of  the  exponential 
dependence  of  pulse  loss  thickness  upon  h*  when  it  is  being  calculated  it  is  possible 
to  disregard  the  second  ecmpoaent  in  comparison  with  the  first  one*  As  was  shewn  by 
calculations*  at  greater  values  h  the  err  or  ^obtainable  when  calculating  the  friction 
coefficient*  is  insignificant  at  such  an  approximat ion* 

Having  accepted  for  pulse  loss  thlckneso  such  an  approximation*  ve  substitute  its 
expression  in  the  Integral  formula  of  pulses.  (1*1)*  After  simple  transform  ve  will 
obtain  the  following  equations 


where 


Taking  the  integral  from  both  parte  of  this  equation  and  utilizing  the  boundary 
condition  Rg  ■  0  at  h  •  0,  va  will  have  in  this  approximation 

xa 

Q.il) 

Log&r itbming  (1,21)  and  turning  away  fraa  variable  h  to  of  according  to  equation 
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we  will  obtain  the  known  Karman  f  oral  la 


{1.3.3) 

The  Just  now  explained  order  of  deriving  the  harman  r  or  mu  in  is  not  in  itself  itnpor* 
tent*  but  only  ea  en  explanation  to  a  more  complex  deriTstloa  of  Karmen  fonxile  genersU 
izatioa  for  the  ease  of  a  gas  flow* 

2#  Turbulent  boundary  layer  in  the  flow  of  **  compressible  gkm 
We  will  confine  ourselves  to  the  case  when  the  Prendtl  number  equals  one  (l*«l) 
and  will  discuss  the  flow  around  a  plate,  as  thermally  insulated  ("lamellar  thermo* 
meter")  and  with  discharge  of  heat  as  veil* 

We  will  assume*  that  in  any  given  section  of  the  boundary  layer  is  valid  the  fol- 
lowing  known  ratio  between  the  absolute  temperature  of  the  gaa  T  and  its  longitudinal 

velocity  m 


where 


n_ , 

<»>  *  1  •  0  when  the  surface  la  heatad 

Jsrl  9  3a — 

and  Uj  /m  0  whan  cooled)  |  genua  *  2  Ki»  Tw  | 


Tt  ■  T/  (1  ♦  2  >£•#  )  -  temperature  of  the  thermometer,  equalling  tba  temper¬ 

ature  of  plate  surfaoe  TW  la  the  absence  of  beat  transfer  («->}■  0). 

Noting,  that  at  a  praeaur#  constancy  ewer  tba  profile  of  the  layer  takea  place  an 


equality 


i-*-— (*M*f  &x) 


wa  will  eccapile  analogous  with  (1*8)  the  baale  Kerman  fornula,  which  In  tba  presenea 
of  ecopreosibility  of  the  medium  will  have  the  form  ofi 

Substituting,  oo  boforo  in  par,l«  tha  arguaant  and  fu notion  plseos,  wa  will  obtain 
instead  of  (2*3)  an  agnation 

as  first  integral  in  whiah  la  a erred  by 


YyJ* - 

~  I  Xh  .  r  r  h  I>  J  Y 

r]  -  C  v\p  I  l:c  sin - J- 

V  Y  _ ’ 

<V 


1 


(M) 


To  determine  the  cojm tents  C  it  is  necessary  to  deal  with  the  value  of  the  derive 

!£. 

tive  on  the  boundary  of  the  lamlnary  sublayer  from  the  side  of  the  turbulent 

nuoleus,  A  simpler  assumption  appears  to  be  the  requirement  that  this  value  should 

_L. 

be  of  the  very  same  at gnitude  f  •  xa  ,  as  in  the  ease  of  low  velocities}  La  other 
words,  it  is  assumed,  thaty  (a*0)  does  not  depend  upon  the  guana  and  omega  coefficients 
Such  an  assumption  is  equivalent  to  disregarding  a  change  in  density  along  the  length 
of  the  lamlnary  sublayer;  it  leads  to  the  following  results 


il 


]  \ 


i. rc  b;n - 

,v  \  i 


_a«» 


]  \  ‘  4- 


/y 


M 


One  can  be  convinced  easily,  that  at  omega  *  garana*  0  the  right  part  of  equation 
(2.6)  coincides  with  the  corresponding  aquation  (1.12)  for  incompressible  fluid# 

After  this  we  will  take  up  the  thickness  or  pulse  loss  of  the  boundary  layer.  In 
case  of  a  compressible  gas  this  value  is  expressed  by  equations 


t(n) 


Using  aquations  (2,2)  and  (2,6)  and  making  transforms,  analogous  as  in  previous 
paragraph,  we  will  obtain  the  following  expression  for  the  thiekness  of  pulse  loss* 


(1  —  U)  —  v)  h'\^ 

/bee 


,  f  «n  —  u) 
j  1  —  uu  —  Y"J 


x  exp 


-4=  I  arc  sin  - 


Vv «  - 


2  V  Y 


V\  a-r 


1 - arc  sin- 


4V 


Introducing  •  new  variable 


we  transform  (2*8)  into  fora  of 


jrv 

Vv  “  *r 

^  =  tt=  arc  sin - T- 

V Y  1/ . 


(*,«> 
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where 


<1—  <u  —  y>  v.*»  !•  ] 


r — "t  /* 


Ji_  )  '2L>*±y 

V  ^I'W 


x  ( bin |/ Y  ^  —  sin l ' y  v,  _  »ir. |  v 't -  ~  sin  |  y  ^'jy. 

y  v  | 

v'  lx'-w 


—  -  ;:rc  bin  — —  ; 

V  V  V  1  H-  -7- 

1  «  -»v 

Vyi--± 


ki.M) 


(2.1  o) 


Yj  '  77  aa*  7— . 

VY  V'*%  j 

To  calculate  the  determined  integral,  situated  in  the  right  eide  of 
equation  (2.10),  we  will  again  assume  that  the  method  of  symtotio  analysis, 
the  same  as  in  the  previous  paragraph.  Intergrating  with  an  accuracy  to  the 
member  j  containing  x^h^  in  the  denominator  and  substituting  approximation^^ 
in  form  of  V,  -  V.  £  JX 


*-  -  v,  -  -  v  • 


we  will  obtain  the  following  expression  for  pulse  thickness  loss 

s»*  c~*aii-  **  - Y)  v.  S'- ^ i' i — ?rJ5,"_v  ( a  | a) 

6  ***<t  x/<  n-«-vi  \X  J  7/ 

In  an  analogous  way  it  Is  possible  to  obtain  the  expression  for  the 

i  -  1 

thickness  of  displacement  Ylv  i  y  (3- y - --)  I 


we  will  find i 


thickness  of  displacement  _  ,-x« (I  .y,Ve  ,  _ 21— i -v  1  I  •  Iq  ,u\ 

°  "  /x-t/.  L  wm -i-r* ->*-y  j  \i*'V 

Compiling,  according  to  equations  (2.14)  and  (2.13),  the  ratio  c)  */j>**«  H, 

we  will  findi  y(3~  I-5u>— Y  — 'o^r) 

_  **  ~~ v)  | 1  ~  y.ft  (i  —  y>  j  (j  \s) 

_  (l_“"”Y)l  xA/l-w-y] 

In  a  epee if io  case  of  flow  of  an  incompressible  fluid  (gammaaO)  and  in  the 
absence  of  heat  transfer  (omega  •  o)  we  obtain  a  term  for  1,  perfectly 
coinciding  with  the  corresponding  equation  of  the  previous  paragraph  (1.19)* 

Having  assumsd,  as  in  previous  paragraph,  for  pluse  loss  thickness  a  more 
cruder  approximation  (disregarding  the  seeond  component  in  eq.  (2.13)  we  will 


substitute  this 
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equation 


x  •  *  __  £  Kgv^(l-M  — Y>  e»MV.-»») 

into  an  intagral  aquation  of  pulaesf  coinciding  in  fora  with  (1#1)#  After  simple  trans. 
fonna  we  will  obtain  an  equation* 

(X.n) 


where 


Having  taken  the  integral  from  both  parts  of  Jthis  equation  and  using  boundary 
condition  *  Oath*  0f  we  will  obtain  in  this  approximation 

-C r  dV*  =  Rss  (p,  I %) 

Equation  (2#l8)  transforms  at  the  limit  into  a  corresponding  expression  for 
inccmpreasible  fluid  (g  fi) 

Next  we  vill  separata  both  parts  of  equation  (2.18)  into  eonfornal  two  parts  la 
equation  (219  )•  and  then  we  will  have 

\irj  e  o-t‘-  vJL.JUv 

Changing  over  in  (2«20)  from  h  to  Of  and  logarithming,  we  will  obtains 


_iciZ!L 

cl  el» 


('(’«  —  tf  J  V 1  —  w  —  V 


/i 


-“'fJ-  (2.. 2,1) 


Accepting  for  Of0  a  gradual  dependence  upon  the  number  Rj* 

c,„  =  0,0263/?- v* 

7  *0 

and  introducing  designations 

f  =  „  l,55/?ioM; 

Vc/> 


15,. n) 


1  /]  —  to  —  Y 

f  ^ 

arc  c?n  _ 

»  Y 

fli  t  dill 

1  1 

2  Ky 


_ _ = - arc  sin  - — 7  , . 


fa) 

Wy 


ws  will  obtain  s  dependence  of  the  frietion  coefficient  upon  the  parameters  of  ths 
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unocmlng  flow  and  condition*  on  the  wall  in  the  following  unclear  forms 

(a-2+) 


Equation  (2.24)  can.  be  reduced  to  a  tranacedental  equation  with  one  parameter  In 

form  at  i 

i  e.v-f.v-L,  (J.2  5j 


where 


*n  this  way,  the  expression  fer  the  friction  coeff  icientjcan  be  vrittens 

^  _  fzKz 

cu  ~  ‘ 

where  F  an  K  -  known  functions,  end  N  is  determined  from  solutic  n  of  equation  (2*25). 

As  was  shown  by  calculationaf  equation  (2,25)  can  be  replaced  by  an  approximated 
one,  by  a  simpler  equations 


.V  -  0,123  -r  0,820  L, 


(JL  %%) 


which  appears  to  be  sufficiently  accurate  in  a  wide  range  of  changes  of  all  paramo » 

ters.  In  this  case  equation  (2*27)  acquires  the  form  ofs 

c ,  r-A'* 


<  [o.ias-, 0.620 (!g^  S-r-))  (7L&$ 


To  determine  value  G  it  is  necessary  to  accept  a  law  of  change  in  viscosity  due 
to  temperature*,  In  the  role  of  such  a  law  can  be  accepted  either  the  Sezerlend  formic 
la,  or  the  gradual  viscosity  dependence  upon  temperature*  In  the  latter  case  the  ex* 
pression  for  G  has  ths  form  of  t 


.«->*(%)•  (2.30) 

where  a-  index  of  the  d  grea  in  thaj^/-'la  law.  It  would  be  interesting  to  compare  the 
described  method  with  the  Van-Drayat  method,  which,  as  is  known,  is  bassd  on  ths  trans¬ 
fer  of  the  Brand tl  formla  into  gaa  dynamics.  After  the  Van-Hrayst  formula  is  brought 
to  form  (2.24)  it  beeoass  explained,  that  the  difference  lies  in  the  Talue  0.  which 
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during  th#  adoption  of  the  gradual  law  of  change 


in  viscosity  is  expressed  by  equa* 

l.iOoj  . 


ri(.u  ?i«*2 

Upon  ths  examination  of  equation  (2.29)  it  is  revealed  that  at  greater  values 
of  Re-number  the  magnitude  of  the  ratio  Cf/efo  depends  very  weakly  upon  the  Re-number. 
The  truth  is.  directing  the  Re-number  into  infinity  (Rj0  )  or  I  toward 

infinity,  we  obtain  the  following  specific  formula  for  the  ratio  Of/e^t 

£-**•  0 *■) 

This  formula  may  appear  to  be  useful  when  making  approximate  calculations  of 
friction  resistanoa  at  greater  Re- numbers. 

According  to  the  described  method  (fornula  (2.29((  the  local  friction  coefficient 
was  calculated.  Ifce  results  of  calculation  for  the  ease  of  absence  of  heat  transfer 
(1^/T^  ■  1)  and  ccoperieon  of  same  with  experimental  data  of  different  authors^* 

|8]  are  shown  in  fig.l^Three  curves  on  this  drawing  correspond  to  three  different  values 
of  Re-number  (lO^.lO^.lO®)* 
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As  la  evident  from  tha  drawing,  the  calculated  and  experimental  data  are  la 
excellent  conformity  all  the  way  up  to  greater  M  number  a. 

The  effect  of  the  temperature  factor  on  the  friction  coefficient  la  shown  in  fig.2. 
We  shall  examine  two  available  practical  values  of  partial  cases* 

1)  Plate  Heat  Insulated 

In  this  ease  the  temperature  of  the  wall  is  equal  to  the  temperature  of  de— 
celeratioat 

t.  M1  -T-  '1--)-  QUH) 

and  the  values  anega  and  gaum*  respectively  are  equal  to 


Function  K  determining  the  process  will  bet 

1  arc  sin  1  7.  (^,33 ) 

2)  rlate  with  heating  or  cooling  at  small  velocities 
In  this  case9  guiding  gamma  we  will  haves 

K..JS1CIEF-- ■>.  (2_3fJ 

In  the  special  case  of  flow  of  an  incompressible  fluid  in  the  absence  of  heat 
transfer  (  omega  *  ganma  =  0)  the  value  K  *  1* 
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